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We present a sufficient condition under which a weak solution of the Euler- 
Lagrange equations in nonlinear elasticity is already a global minimizer of 
the corresponding elastic energy functional. This criterion is applicable to 
energies W(F) = W(F'^F) = W (C) which are convex with respect to the 
right Cauchy-Green tensor C = F'^F, where F denotes the gradient of de¬ 
formation. Examples of such energies exhibiting a blow up for det E —)• 0 are 
given. 
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1 Introduction 

In a recent contribution, we have stated conditions for an elastic energy W{F) = W{C) 
to be convex in terms of the right Cauchy-Green deformation tensor C = F^F. These 
conditions include dependencies on the logarithm of the determinant, showing the nec¬ 
essary singular behavior for det E ^ 0; see m for more details. 

This subject has been taken up by Spector m and Silhavy [16] , who presented alter¬ 
native proofs. The main result is the following. Here and throughout, Sym+(n) denotes 
the set of positive definite symmetric n x n-matrices and M+ := (0, oo). 


'^School of Applied Science, Federation University Australia, Mt Helen VIC 3350, Australia, email: 
d.gao@federation.edu.au 

^Patrizio Neff, Head of Lehrstuhl fur Nichtlineare Analysis und Modellierung, Fakultat fiir 
Mathematik, Universitat Duisburg-Essen, Thea-Leymann Str. 9, 45127 Essen, Germany, email: 
patrizio.neff@uni-due.de 

^Department of Mathematics, University of Craiova, A. I. Cuza Street, Nr. 13, 200585 Craiova, Roma¬ 
nia, email: ionelroventa@yahoo.com 

■^Corresponding author: Christian Thiel, Lehrstuhl fur Nichtlineare Analysis und Modellierung, 
Fakultat fur Mathematik, Universitat Duisburg-Essen, Thea-Leymann Str. 9, 45127 Essen, Ger¬ 
many, email: christian.thiel@uni-due.de 


1 



Theorem 1.1. Let f G C'^(]R+). The function C i-)- /(detC) is convex in C on the set 
Sym+(n) if and only if 

71 — 1 

f''{s) H- f'{s) > 0 and f'{s) < 0 for all s G M+ . (1) 

n s 

In particular, these conditions are satisfied for f{s) = — logs. 

In this short note we would like to substantiate our claim from m that convexity 
in C should be somehow “nice”. It is clear that convexity of the strain energy in C 
is independent of both the rank-one convexity condition m and Ball’s polyconvexity 
condition [T]. Therefore, convexity of W as function of C is not sufficient for establishing 
the existence of minimizers via the direct methods of the calculus of variations. Of course, 
in a neighborhood of the identity, an existence proof based on the implicit function 
theorem is always possible [3]. 

However, in m we have given an example of a strain energy W{F) = W{C) which is 
both polyconvex with respect to F and convex in C] furthermore, the energy shows the 
correct behaviour for infinite compression, i.e. W(C) —t- oo as det C —t- 0. Convexity with 
respect to C is also used profitably in the derivation of some models in elasto-plasticity, 
see [8]. 

2 Applications of the convexity with respect to C 

We consider the problem of minimizing the elastic energy I(ip) = Jq Wif^ip) dx, where 
<y9: H ^ is the deformation of an elastic body H C M^, subject to the boundary 
condition </5|r = (po with T C dLl, see [2]. The corresponding formal Euler-Lagrange 
equation is 0 = DivDi7’iy(T) = Div5i(T), where F = 'V(p denotes the gradient of de¬ 
formation, Si{F) = DpW{F) = F S 2 {F'^F) is the first Piola-Kirchhoff stress tensor and 
S 2 iC) = 2 DcIT(C') is the second Piola-Kirchhoff stress tensor. The following simple 
observation extracted from mm shows the applicability of our convexity condition. 

Proposition 2.1. Assume that W is convex in C and that po is a weak solution of 
the corresponding Euler-Lagrange equation. If the second Piola-Kirchhoff stress tensor 
S 2 {Co) = 2 DcIT(Co) corresponding to Cq = is positive definite everywhere in 

LI, then this solution is a global minimizer. 


Remark 2.2. Note that the condition of positive definiteness only refers to the second 
Piola-Kirchhoff stress evaluated at the weak solution pQ of the Euler-Lagrange equation. 

Convexity with respect to C is tantamount to the monotonicity of the second Piola- 
Kirchhoff stress tensor S 2 , i.e. to the condition 

{S2{Ci)-S2{C2), Ci-Cs) > 0, 

where {X,Y) = tifY'^X) denotes the canonical inner product on This does not 

violate any physical principle, whereas monotonicity of the first Piola-Kirchhoff stress 
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tensor 5i (i.e. convexity with respect to the deformation gradient F) must be excluded 
to ensure physically plausible material behaviour. 

The positive definiteness of the second Piola-Kirchhoff stress 52 (Cq) at ipo, on the 
other hand, means that all occurring forces, transformed to the reference configuration, 
are of a “tensional” type, i.e. that the solution (pQ must not exhibit any “compressional” 
forces. Note, again, that this condition only refers to the solution cpQ and that we do not 
require 82 ( 0 ) = 2 DchT(C') to be convex for all C G Sym+(3), which together with the 
convexity of VT in C would imply convexity of W in F. 

Proof. We write (p as (p = ipo + u. Then Vp> = V(/?o + Vu and thus 

Vpy^Vip = (Vy^o + Vu)'^(V(/Jo + Vn) = VipqVp>o + Vcp^Vu + Vu^Vipo + . 

We therefore find Vpf^Vip— 'VpiQVipo = Vp>qVu + Vu'^Vp>o + Vu^Vu. The convexity 
of W with respect to C implies 

WiC) - W(Cq) > {BcW{Co), C - Co) for all C G Sym+(n). (2) 

Using ([2]) we obtain the following minimizing estimates: 


in 


W(V(^^V<y9) - W(V(^JV^o)dx 


> / (Dc'W(V(/?q V(/9o), V(/?Q Vtt + Vu^V(/Jo + Vn^Vu) dx 

Jn 

= [ 2{Si{Vipo),Vu) + l-{S2{Co),Vu^Vu)dx 
Jn ^ 


I [ (52(Co),Vu^Vn)dx > 0, 
^ Jn 


(3) 


where the equality and inequality in ([3]) follow from the fact that 2(Si(V(po),Vu}dx = 
0 due to (po satisfying the weak form of Euler-Lagrange equation and the assumed positive 
dehniteness of 52 ( 6 * 0 ), respectively. ■ 


If the convexity condition in Proposition 12.II is strengthened to the uniform convexity 
requirement DqW{C).{H,H) > c'^ \\H\\‘^ for all symmetric matrices H and some c"*" > 0, 
and if the solntion ipo is a diffeomorphism, then it is easy to show that 

j^DlW{Vipo).{VfJ^Vf))dx > c-^\m\hin) 


for all test functions iJ G C^{Q) and some c+ > 0, implying the local stability of the 
solution [T^ . 

We will apply Proposition 12.II to a modified version of W as given in |11| . 
Proposition 2.3. Let 

W{F) = W{C) = W{F^F) = a{tr{C))^ + /3tr(C2) - logdetC 
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OTt/i 6a+2/3 — 1 = 0 (which is strictly convex in C and poly convex provided that a, (3 > Oj. 
Assume that W{C) admits a smooth solution +o of the Euler-Lagrange equation satisfying 
detV(/?o > C~^ > 0 and S' 2 (V+o) is positive definite everywhere. Then cpQ is a global 
minimizer. 

The condition 6a + 2/3 = 1, a,/3 > 0 ensures that the condition DcVh(l) = 0 
for a stress free reference conhguration is fulfilled. Similarly, for the energy W{C) = 
atr(C') + /3tr(C'^) — logdetC, we need to impose the conditions a + 2/3 = l, a>0 
and /3 > 0. Since the energy W given in Proposition 12.31 is polyconvex and coercive 
in the Sobolev space W^’'^ if a > 0, it follows from the direct methods of the calculus 
of variations [2] that there exists a global minimizer to the boundary value problem 
considered in Proposition 12.11 The conditions under which a global minimizer satishes 
the weak form of Euler-Lagrange equation are, however, not clear since it is not known 
whether det V+o ^ c"*" for some > 0. On the other hand, most computational methods 
will try to solve the Euler-Lagrange equations. Thus, in general, it is not known whether 
a global minimizer has been found. 

The energy 

W{C) = /i(a(tr(C))^ + /3tr(C) -logdetc) 

is a special member of the compressible Neo-Hooke family. It generates the linear elastic 
response p, || devep + /r ^2 a + |^ [tr(e)]^ with an induced linear Poisson coefficient u = 

h ia+i ^ (Oj ^)- Here, || . || denotes the Frobenius matrix norm and e is the linearized 
strain tensor. Our examples show that the search for further insight into the properties 
of specihc elastic energies might lead to unexpected new results. Certainly, the class 
of convex energies in C is much to narrow from an application point of view. Since, 
however, all invariance principles and side conditions have been respected, our examples 
show that there may be enough space for more discoveries. 

A more general application of the convexity in C was given by Le Dret and Raoult m- 
they showed that if W{F) = W{C) is convex with respect to C, then the quasiconvex 
hull QW of W is given by an analogue of Pipkin’s formula |13j : 

QW(F) = inf W(F^F + S) . (4) 

SeSym+(3) 

3 Further examples 

A trivial example of a function W which is convex in C is the Saint Venant-Kirchhoff 
energy MUM 

^Hsvk(C’) = + ^tr(C-l)^ 

which is convex in C as long as /i > 0, 3A + 2/i > 0. However, I+svk is not rank-one 
convex, as has been shown by Raoult |14j . 
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For another class of examples, consider an energy function satisfying the Valanis- 
Landel hypothesis, i.e. an energy of the form 

W (F) = w{Xi) + w{\ 2 ) + w{X 3 ) 

where w: M+ —)• M is a fixed function and Ai,A 2 ,A 3 are the singular values of F, i.e. 
the eigenvalues of U = \fC. Since the eigenvalues of \J are the square roots of the 
eigenvalues of C, the function W can be represented in terms of the eigenvalues of C via 
W{F) = w{Xi)+ w{X 2 )+ w{X‘^) with w{t) = w{y/t). Then according to a result by Davis 

[3], VF is convex with respect to C if and only if w is convex. If /; M_|_ —>■ M satisfies the 
requirements of Theorem 11.11 then an energy of the form 

1F(F) = u;(Af) +'«;(A|) + ^(AI) + 5f(detC') 

is convex in C as well if w is convex. 
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